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ABSTRACT 

We present the first f ar-IR observations of the solar- type stars 6 Pav, HR 8501,51 Peg and f 2 Ret, taken within the context of the DUNES Herschel 
Open Time Key Programme (OTKP). This project uses the PACS and SPIRE instruments with the objective of studying infrared excesses due 
to exo-Kuiper belts around nearby solar-type stars. The observed 100 yum fluxes from 6 Pav, HR 8501, and 51 Peg agree with the predicted 
photospheric fluxes, excluding debris disks brighter than L& mt /L* ~ 5 x 10~ 7 (1 cr level) around those stars. A flattened, disk-like structure with a 
semi-major axis of ~ 100 AU in size is detected around £ 2 Ret. The resolved structure suggests the presence of an eccentric dust ring, which we 
interpret as an exo-Kuiper belt with Ld llst /L* ~ 10~ 5 . 
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1. Introduction 
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The discovery of infrared excess emission produced by cold, op- 
tically thin disks composed of micron-sized dust grains around 
main sequence stars is one of the main contributions of IRAS 
(Aumann et al. I1984I I. Since the lifetime of such grains, set 
by destructive collisions, Poynting-Robertson drag and radiation 
pressure, is much shorter than the ages of these stars, one must 
conclude that those dust disks - called debris disks - are continu- 
ously replenished by collisions of large rocky bodies (Backman 
& Paresce [1993l ). Observations of debris disks provide power- 
ful diagnostics from which to learn about the dust content, its 
properties and its spatial distribution; in addition, since dust sen- 
sitively responds to the gravity of planets, it can be used as a 
tracer of the presence of planets. Thus, observations of debris 
disks around stars of different masses and ages inform us about 
the formation and evolution of planetary systems, since they are 
a direct proof for the existence of planetesimals and an indirect 
tracer of the presence of planets around stars. 

In the Solar System, the asteroid and Kuiper belts are exam- 
ples of debris systems; in particular, the Kuiper belt has an esti- 
mated dust luminosity L dust /L* ~ 10' 7 -10~ 6 (Stern [T99oT l. IRAS 
was only able to detect bright disks, Ldust/^* > 10~ 4 , mainly 
around A and F stars; ISO extended our knowledge to a wider 
spectral type range and found a general decline with the stellar 
age (Habing et al 2001 Decin et al. 2003). A remarkable step 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



forward has been achieved by Spitzer, studying debris disks as 
faint as L,j ust /L* several times 10~ 6 , their incidence from A to M 
type stars, the age distribution, the presence of planets, etc. (e.g. 
Su et al. 120061 Trilling et al. 2008] Bryden et al. 120091 . Spitzer 
has, however, several limitations. Its poor spatial resolution pre- 
vent us from constraining fundamental disk parameters which 
require resolved imaging, and the confusion limit inherent in its 
large beam limits its detection capability to cold disks brighter 
than the Kuiper belt by around two orders of magnitude. Also, 
Spitzer is not sensitive longward of 70 pm, wavelengths partic- 
ularly important for the cold disks generally found around Sun- 
like stars. The far-infrared 3.5 m diameter Herschel space tele- 
scope (Pilbratt et al. 2010) with its instruments PACS (Poglistsch 
et al. 2010) and SPIRE (Griffin et al. 120101 ) overcomes these lim- 
itations, offering the possibility of characterising cold, ~ 30 K, 
debris disks as faint as Ldust/^* few times 10~ 7 with spatial res- 
olution ~ 30 AU at 10 pc, i.e., true extra-solar Kuiper belts. 

DUNESQ is a Herschel OTKP designed to detect and charac- 
terize cold, faint, debris disks, i.e., extra-solar analogues to the 
Kuiper belt, around a statistically significant sample of main- 
sequence FGK nearby stars, taking advantage of the unique ca- 
pabilities of Herschel with PACS and SPIRE. The data will 
be analysed with radiative, collisional and dynamical dust disk 
models. A complete description of DUNES goals, target selec- 
tion, and stellar properties will be presented elsewhere (Eiroa et 
al. 2010, in preparation). The objectives of the DUNES survey 
are complementary to those of the OTKP DEBRIS (Matthews et 
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Table 1. Summary of the SDP DUNES observations. 



Star 



Obs. ID 



Mode 



Bands 
(A™) 



Scan 
Direction 



OT 

(sec) 



q 1 En 
q 1 Eri 
q 1 Eri 
<T 2 Ret 
<f 2 Ret 
<T 2 Ret 
(5 Pav 
HR8501 
51 Peg 



1342187139/40 

1342187141 

1342187142 

1342183660 

1342191102/03 

1342191104/05 

1342187075/06 

1342187145/46 

1342187255 



SM 

PS 

PS 

PS 

SM 2 

SM 2 

SM 

SM 

PS 



100/160 
100/160 

70/160 
100/160 
100/160 

70/160 
100/160 
100/160 
100/160 



637117° 



117763° 
63717° 
457135° 
637117° 



2536 
4714 
789 
1572 
4510 
4510 
3834 
1890 
1731 



(1) SM = scan map; PS = chop-nod/point-: 

(2) Routine phase DUNES observing time 



source mode 
(not SDP) 



al. 120101 . Both projects have complementary star samples, shar- 
ing partly some sources and the corresponding data. 

The DUNES objectives require the detection of very faint 
excesses at the mJy level, comparable to the photospheric emis- 
sion and only a few times the measurement uncertainties. The 
primary observing strategy is designed to integrate for as long as 
needed to detect the 100 //m photospheric flux, subject only to 
confusion noise limitations. In this letter we present our first re- 
sults obtained during the science demonstration phase (SDP) of 
Herschel. Four solar-type G stars were observed: £ 2 Ret, 8 Pav, 
HR 8501, and 51 Peg. We also observed q 1 Eri as a test object 
with a well-known, bright debris disk; the q 1 Eri results are pre- 
sented in an accompanying letter (Liseau et al. 1201 Oi l. Excluding 
8 Pav, the rest of the stars are shared targets with DEBRIS. 

2. Observations and data reduction 

The stars were observed with PACS at 70 /mi (blue), 100 /mi 
(green), and 160 /mi (red). Two observing modes were used - 
chop-nod/point-source (PS hereafter) and scan map (SM here- 
after). Our data set allows us to make a direct comparison of 
both modes, specially in the cases of q 1 Eri and ( 2 Ret. A criti- 
cal evaluation of this comparison will be the subject of a techni- 
cal note. SM observations of £ 2 Ret were carried out as DUNES 
routine observations. Table 1 gives some details of the observa- 
tions including the identification number, the observing mode, 
the wavelength bands, the scan direction angles, and the total 
duration of the observation (OT). 

Data reduction was carried out using the Herschel Interactive 
Processing Environment (HIPE), version v2.0.0 RC3, and the 
pipeline script delivered at the December 2009 Herschel data 
reduction workshop, ESAC, Madrid, Spain. The script provides 
all the tools to convert pure raw PACS/Herschel data to flux cal- 
ibrated and position rectified images. While the instrument pixel 
size is 3'.'2 for the blue and green bands and 6'.'4 for the red band, 
the resolution of the final images is set to 1 "/pixel and 2"/pixel 
for the blue/green and red bands, respectively. 

3. Results 

3.1. 8 Pav, HR 8501 , 51 Peg 

Table 2 gives J2000.0 equatorial coordinates of 8 Pav, HR 8501 
and 51 Peg at 100 /mi, as well as their optical positions. PACS 
positions are corrected from the proper motions of the stars. 
Differences between the optical and 100 /mi positions are within 
the uncertainties for Herschel pointing. Of the three stars, only 8 
Pav has been detected at 160 /mi. 



100 /mi FWHM values of 8 Pav have been estimated us- 
_ing a 2-D gaussian fit. This procedure did not produce rea- 
sonable results for HR 8501, perhaps due to the faintness of 
_the star; in this case the FWHM has been estimated from in- 
tensity profiles along the R.A. and Dec. directions. For 51 
Peg, observed in PS mode, the FWHM estimate is also based 
on some additional point-like sources visible in the PACS 
field. The 100 /mi FWHM estimates are given in Table 2. 
The 160 /mi 2-D gaussian fit for 8 Pav gives FWHM = 
ll"8x9"2, with conservative errors ~ 1". These values and 
the elongated beam are consistent with the expected results for 
point sources (see the technical notes PICC-ME-TN-035/036 in 
"http://herschel.esac.esa.int/AOTsReleaseStatus.shtml) 

Aperture photometry has been used to estimate the flux of the 
stars. Table 2 gives fluxes and errors taking into account the cor- 
rection factors indicated in the aforementioned technical notes. 
The sky noise was 2.5xl0~ 5 Jy/pixel and 2.7xl0~ 5 Jy/pixel for 
the 100 /mi SM data of 8 Pav and HR 8501, respectively. The 
sky noise is considerably higher in the PS-mode 51 Peg image 
(w4.3 x 10~ 5 Jy/arcsec 2 ) due to the very irregular background 
and presence of negative signals. The 160 /mi sky noise was 
4.9x10 Jy/pixel for the 8 Pav data. The absolute calibration 
uncertainties are 10% in the blue and green bands and better 
than 20% in the red band. 

3.2. ( 2 Ret (HIP 15371) 

Fig 1 shows the SM PACS images of £ 2 Ret. An East- West ori- 
ented structure is seen at 70 and 100 /mi. It consists of two point- 
like flux peaks embedded in a faint, extended emission, which 
displays a secondary diffuse maximum at its Western side. Both 
point-like peaks have similar brightness in the green band, but 
the Eastern point-like peak is much fainter in the blue band. The 
two point-like sources are unresolved in the lower-resolution 160 
//m image; instead a single bright peak is observed at that posi- 
tion with a secondary maximum at the position of the 70/100 /mi 
Western diffuse emission. 

Table 3 gives positions at 70 and 100 //m of both point-like 
sources, and of the brightest 160 //m peak; the optical position 
of £ 2 Ret is also given for comparison. The brightest 70 /mi peak 
coincides with the optical position of the star within the Herschel 
pointing error; this result and the fact that its PACS 70 and 100 
//m fluxes are similar to the expected ( 2 Ret photospheric fluxes 
(see below) lead us to identify this point-like PACS object with 
the optical star. There is a small shift between the 70 and 100 /mi 
positions of f 2 Ret, but we note that a similar shift is found for 
other field objects - a blue object very close to the ( 2 Ret com- 
plex towards the South- West; and two red objects, one towards 
the Northwest and one towards the Northeast (see Fig. 1). 

The middle panels of Fig. 1 show isocontour plots. The op- 
tical position of f 2 Ret is marked. 100 //m and 160 /mi contours 
have been spatially shifted so that the peak positions of the men- 
tioned field objects coincide in all three bands (those objects are 
not shown in the isocontour plots). The size of the whole struc- 
ture changes with wavelength from » 25" x 15" in the blue to 
» 40" x 15" in the red band. East-West intensity profiles are 
shown in the bottom of Fig. 1 , together with similar profiles of a 
Bootis. The blue and green intensity profiles show the point-like 
character of Z 2 Ret, as well as of the faint peak at the East, called 
PS-E hereafter ; the profile in the red band also shows point-like 
behaviour for the brightest 160 /mi peak, £ 2 Ret+PS-E in Fig. 1 
and Table 3. The Western diffuse peak ("W" in Fig. 1) appears 
very prominent in the green and red profiles, while it is very faint 
compared to f 2 Ret in the blue profile. North-South profiles (not 



C. Eiroa et al.: Cold dust around nearby stars (DUNES). First results 

Table 2. Equatorial coordinates, FWHM at 100 yum, observed fluxes with lcr statistical errors (Fpacs), and predicted photospheric 
fluxes (F+). Absolute PACS uncertainties are ~10% and less than 20% for 100 and 160 /urn, respectively. Flux units are mJy. 



Star 



Optical position 
(J2000) 



PACS 100 yum position 
(J2000) 



PACS 100 fim 
FWHM Fp A cs 



PACS 160yum 

fpACS F* 



cSPav 20:08:43.61-66:10:55.4 

HR8501 22:18:15.62-53:37:37.5 
51 Peg 22:57:27.98 +20:46:07.8 



20:08:43.53-66:10:58.1 ft' 3x67 2 59.6+1.1 68.7 
22:18:16.14-53:37:37.1 6'.'3x6'.'2 9.8+1.2 10.9 
22:57:28.44+20:46:08.5 5'.'9x6'.'0 11.3+1.7 10.8 



21.0+1.5 



26.9 
4.3 

4.2 



70p,m 



100 urn 



PS-E I W 
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Fig. 1. PACS results of £ 2 Ret. Panels from left to right: 70 yum, 100 yt/m, and 160 /jm. Upper panel: Field size is 100" x 100" with 
East to the left and North up. Middle panels plot isocontours. Note that the field size is different from that displayed in the upper 
panels. The "star" symbol, position (0,0), corresponds to the optical star (Table 3). A segment indicating 120 AU is shown. Contours: 
70yum: 3,4,6,9,12,15,24 cr, 100 yum: 3,4,5,6,12,14 cr, 160 fim: 3,6,9,12,15,18 cr ( lcr values are indicated in the text). Lower panels: 
R.A. intensity profiles through the peak flux of ( 2 Ret (solid line). For comparison, R.A. intensity profiles of the calibration star a 
Bootis are also plotted (dashed lines), scaled to the peak flux of £ 2 Ret (70 and 100 fim) and I 1 Ret+PS-E (160 fim). The a Bootis 
profile is also scaled and superimposed to PS-E in the 100 fim panel. 



shown) do not resolved either £ 2 Ret, PS-E, or £ 2 Ret+PS-E in 
any band, i.e. they are point-like along that direction with no hint 
of any faint extended emission. 

Table 3 gives PACS fluxes estimated using the flux peaks 
of the point-like sources and integrating over beam sizes given 
by n (FWHM A - x FWHM,,) /4 In 2. The PACS source identified 
with ( 2 Ret is a blue object, while PS-E is a red one. The flux at 
160 fim corresponds to ( 2 Ret+PS-E, but PS-E is the main con- 
tributor to the flux at this wavelength - the emission is peaking 
more towards PS-E (Fig. 1). The total flux of the £ 2 Ret com- 
plex is 44.5 mJy, 40.4 mJy, and 42.6 mJy in the blue, green and 
red bands, respectively. The estimated 70 yum flux for the whole 



complex agrees very well with the unresolved Spitzer flux of 46 
mJy at the same wavelength (Trilling et al. 2008). 



4. Discussion 

Our data do not reveal any cold dust disk around 8 Pav, HR 8501 
or 5 1 Peg, since the observed and predicted 1 00 fim photospheric 
fluxes coincide within the uncertainties (Table 2). Assuming dust 
temperatures of 40 K (peak blackbody fluxes at ~ 100 fim), we 
can exclude debris disks with Ljust/^* >5x 10~ 7 (lcr) (Eqn. 4 
from Beichmann et al. 2006). 

f 2 Ret, located at 12 pc, is a Gl V star with a bolometric 
luminosity of 0.97 L and an estimated age of ~ 2-3 Gyr (Eiroa 
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Table 3. Coordinates,fluxes, and lrx statistical errors of the ( 2 
Ret PACS point-like sources. Flux units are mJy. 



Object 


<r(2000.0) 


5(2000.0) 


Flux 


( l Ret (optical) 


03:18:12.82 


-62:30:22.9 




C, L Ret (70pm) 
<T 2 Ret (100 pm) 


03:18:12.77 
03:18:13.12 


-62:30:23.2 
-62:30:24.4 


24.9±0.8 
13.4+1.0 


PS-E (70 pm) 
PS-E (100 pm) 


03:18:13.63 
03:18:13.55 


-62:30:24.2 
-62:30:24.4 


8.9±0.8 
13.5+1.0 


< 2 Ret+PS-E (160 pm) 


03:18:13.37 


-62:30:23.2 


19.4+1.5 



10" 



£ iooo 



3 too 



E 1 ' 


1 1 E 

■» f Ret : 


bf 


\ x ps-E - 




\- f 2 Ret + PS-E 


D 


\^ □ Total 
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N. a n a = 


r 
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10 
A (M m ) 
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Fig. 2. SED of ( 2 Ret. Optical, 2MASS, IRAS, and Spitzer 
fluxes are indicated by black symbols. Blue triangles are £ 2 Ret; 
red crosses are PS-E; green triangle is ( 2 Ret+PS-E; magenta 
squares are total fluxes from the £ 2 Ret complex. Error bars are 
smaller than the size of the symbols. The solid line is the best^ 2 
photospheric fit (T eff = 5850 K, logg = 4.5, and [Fe/H] = -0.23, 
which are mean values found using the DUNES discovery tool, 
http://sdc.laeff.inta.es/dunes). The dashed line is a 40 K black 
body normalized at the PS-E 100 pm flux. The deduced 160 pm 
flux from PS-E is also plotted with a red cross (see text). 



et al. 2010, in prep.). Figure 2 shows the stellar SED as well as 
PACS fluxes from PS-E and the whole complex; a PHOENIX 
stellar photosphere (Brott & Hauschildt 2005) is also plotted. 
The agreement between the observed 70 and 100 pm fluxes from 
the blue bright point-source and those predicted by the photo- 
spheric fit, 24.7 mJy and 12.1 mJy at 70 and 100 pm, respec- 
tively, is excellent. This photometry and its positional alignment 
with the stellar position support our claim that the PACS blue 
point-like object is indeed £ 2 Ret. On the other hand, the na- 
ture of the extended structure is intriguing. While coinciden- 
tal alignments with background objects are common in IRAS 
all-sky images, the much higher resolution of Herschel makes 
such juxtapositions unlikely within a targeted survey. Based on 
Spitzer source counts of background galaxies at 70 pm (Dole et 
al. 120041 1. we find that the probability of chance alignment with 
a >20 mJy source within 10" is just 10~ 3 . 

The source PS-E is a red object with a black body tempera- 
ture T(70- 100 pm) * 40 K. We have pointed out that both PS-E 
and C 2 Ret are not resolved at 160 pm and that the flux peak at 
this wavelength is closer to PS-E. If we subtract from the mea- 
sured 160 pm flux the stellar flux predicted for Z 2 Ret (4.7 mJy) 
and make the plausible assumption that the residual flux (14.7 
mJy) originates in PS-E, this 160 pm flux for PS-E is again con- 
sistent with a ~ 40 K blackbody (see Fig. 2). PS-E is clearly 
not stellar; we suggest that it is instead orbiting circumstellar 
dust. The contribution of the extended emission to the total flux 



can be estimated subtracting the point-like sources from the to- 
tal flux reported above. The residual flux mainly corresponds 
to the Western diffuse emission since the point-like sources are 
not resolved along the North-South direction. In this case, the 
remaining flux corresponds to black body temperatures in the 
range ~30-40 K, and the total fractional luminosity from the en- 
tire structure surrounding £ 2 is Ld ust /L* ~ 10~ 5 . 

We have the interesting scenario of a Gl V star surrounded 
by optically-thin 30-40 K emission. This is the temperature 
range expected for black body dust grains orbiting at distances 
~100 AU from the star. This is consistent with the projected lin- 
ear distances from £ 2 Ret to PS-E and to the Western diffuse 
emission of ~70 AU and ~120 AU, respectively. The red im- 
age suggests a flattened, disk-like structure with the star located 
asymmetrically along the major axis, while the blue and green 
images suggest it is ring-like given the flux cavity towards the 
West from the star. We interpret the structure in the PACS im- 
ages as a dust ring surrounding £ 2 Ret. We attribute the observed 
East- West asymmetry to a significant disk eccentricity - e » 0.3. 
Similarly, an offset is observed in the Fomalhaut debris disk with 
e « 0.1 (Stapelfeldt et al. 2004). Maintaining a stable eccen- 
tric ring requires an external driving force such as a shepharding 
planet (Wyatt et al. 1999) and in the case of Fomalhaut the pre- 
dicted planet has been been directly imaged (Kalas et al. 2008). 
The disk asymmetry in the £ 2 Ret system may likewise be the 
signature of an unseen planetary companion. While this is an 
exciting possibility, other forces might also produce disk asym- 
metry. For example, interaction with the ISM is probably respon- 
sible for the strong asymmetry observed around HD 61005, since 
its brightness offset is well aligned with the star's space motion 
(Hines et al. 120071 ). A more profound analysis and detailed mod- 
eling of £ 2 Ret and the suggested Kuiper belt is deferred to a 
future paper. 

5. Conclusions 

Our results show the capability of Herschel/PACS to detect and 
resolve cold dust disks with a luminosity close to the solar 
Kuiper belt, which will allows us to deepen our understanding 
of planetary systems, in particular those associated with mature 
stars. Specifically, our data exclude the existence of cold debris 
disks with Ldust/^* >5x 10~ 7 (lcr) around the solar- type stars 
6 Pav, HR 8501 and 51 Peg. On the other hand, the data show 
that ( 2 Ret is a good example where cold disks around nearby 
stars, very much alike the solar Kuiper belt, can be resolved and 
studied in great detail with the Herschel space observatory. 
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